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Abstract 

Botli continuum and emission line flickering arc phenomena directly associated with the mass accretion 
process. In this work we simulate accretion disk Dopplcr maps including the effects of winds and flickering 
flares. Synthetic flickering Doppler maps arc calculated and the effect of the flickering parameters on the 
maps is explored. Jets and winds occur in many astrophysical objects where accretion disks are present. 
Jets are generally absent among the cataclysmic variables (CVs), but there is evidence of mass loss by 
wind in many objects. CVs are ideal objects to study accretion disks and consequently to study the wind 
associated with these disks. We also present simulations of accretion disks including the presence of a 
wind with orbital phase resolution. Synthetic Ha line profiles in the optical region are obtained and their 
corresponding Doppler maps are calculated. The effect of the wind simulation parameters on the wind 
line profiles is also explored. From this study we verified that optically thick lines and/or emission by 
diffuse material into the primary Roche lobe are necessary to generate single peaked line profiles, often 
seen in CVs. The future accounting of these effects is suggested for interpreting Doppler tomography 
reconstructions. 
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1. Introduction 

Cataclysmic variables (CVs) are close binary systems 
composed by a white dwarf (primary) and a red dwarf 
or subgiant star (secondary). The secondary star fills its 
Roche lobe and matter is transferred to the white dwarf. 
Due to the angular momentum of the system (and if the 
magnetic field of the primary is not strong enough) there 
is the formation of an accretion disk around the primary. 

The fiickering (rapid variability) is observed as stochas- 
tic fluctuations in the emitted radiation, with timeseales 
ranging from seconds to tenth of minutes and ampli- 
tudes from cents of magnitude to more than one mag- 
nitude. Flickering is not exclusive to CVs, but it 
is also observed in symbiotic stars (e.g. CH Cyg, 
Mikolajewski, Mikolajewska & Khudiakova 1990), X-ray 
binaries (Reynolds, Callanan & Filippenko 2007, Malzac 
et al. 2003, Baptista, Bortoletto & Harlaftis 2002) and 
pre-main sequence stars (Kenyon et al. 2000, Clarke et al. 
2005). The flickering cannot be related only to accretion 
disks as it is observed in magnetic systems (Thackeray, 
Wesselink & Oosterhoff 1950); an early model to describe 
the flickering in the polar AM Her was proposed by Panek 
(1980). On a wide perspective, the presence of flickering 
seems associated to the process of mass accretion itself. 

The first CV where flickering was observed is UX UMa 
(Linnell 1949). Since then, many studies were made aim- 



ing to locate the flickering source region on many objects. 
The flickering can be originated in many regions of the sys- 
tem. Warner & Nather (1971) observed that the U Gem 
flickering disappeared during the eclipse, and as the inner 
parts of the disk were not occulted during the eclipse, the 
flickering source was associated with the hot spot. Vogt, 
Kreminski & Podersen (1981) observed that OY Car flick- 
ering persisted during the hot spot eclipse, the same be- 
havior was observed on HT Cas by Patterson (1981) and 
Z Cha by Wood et al. (1986). Home, K. & Stiening (1985) 
identifled the flickering in some eclipsing systems as origi- 
nated in the inner part of the accretions disk, Horne et al. 
(1994) also associated the OY Car flickering with this disk 
region. Bruch (1996) associated the flickering source in Z 
Cha with the region near to the white dwarf, but with 
other flickering sources appearing in other photometric 
states. 

As the CVs cannot be spatially resolved, indirect imag- 
ing methods are used. One of these methods is the 
Doppler tomography (Marsh & Horne 1988), where the 
system emissivity is reconstructed from the observed line 
profiles. Diaz (2001) proposed an extension of the Doppler 
tomography, the flickering tomography, where the line 
profile variability is mapped. In this method, we start 
calculating the variance of the observed line profiles into 
phase bins. Then, the instrumental, orbital, and secular 
components are subtracted from such variance, remaining 
only the intrinsic variance component. The instrumental 
component of the variance is calculated as a combination 
of the Poisson noise plus a Gaussian readout noise. The 
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orbital component contribution is considered negligible if 
the data is sampled into small phase bins. The secular 
component is the remaining long term variability of the 
system. The effect of this later component was found to 
be negligible by replacing the long term average profiles 
by profiles computed from individual runs. From these in- 
trinsic variance line profiles Doppler tomograms are calcu- 
lated, mapping the line emission variability in the system. 
This technique was applied to V442 Oph (Diaz 2001), 
where an isolated flickering source was not identified, and 
V3885 Sgr (Ribeiro & Diaz 2007), where the main flick- 
ering source seems to be associated with reprocessed UV 
radiation on the illuminated face of the secondary star. 

The association between accretion disks and jets is ob- 
served in many astrophysical objects, as young stellar ob- 
jects (YSO), low- mass and high-mass X-ray binaries, mi- 
croquasars and active galactic nuclei. Concerning disks 
around white dwarfs, jets are observed in super-soft bina- 
ries (SSS) and some V Sge systems. Among the CVs, jets 
are not detected, except during the evolution of classical 
novae explosions (Retter 2004), where the jet formation is 
related to the additional energy available due to the nu- 
clear reactions on the primary surface. Soker & Lasota 
(2004) explained the absence of jets in CVs using the jet 
launching thermal model (Torbett f 984), originally devel- 
oped to explain jets in YSO. According to these authors, 
the CVs accretion rates are not high enough to produce 
jets. On the other hand, some CVs have mass loss by 
wind, and the mechanism that best describes the wind 
production in these systems is the line opacity (line driven 
winds) (Proga 2005). 

Ultraviolet observations allowed the identification of 
mass loss by wind in CVs. UV rcssonant lines (e.g. CIV 
A1548, A1551, NV AI239, A1243, Si A1397) arc observed 
with a component in absorption displaced to the blue in 
classical novae remnants and nova-like systems with high 
accretion rate. The fact that wind is detected only in 
systems with high accretion rates relates the wind emis- 
sion to the accretion process. The P-Cygni profiles are 
observed only in systems where the orbital inclination is 
approximately lower than 65.° (eg. Warner 1995; Shlosman 
& Vitello 1993), indicating that the geometry of the line 
formation region could be bipolar. Models from Drew 
(1987) indicate that, even if the emission of the accretion 
disk is bipolar, the wind also must be a bipolar geometry. 

Most of the studies about winds in CVs are based on 
UV ressonant lines. They basically follow two research 
lines: cinematic models fitted to observed profiles (Long 
& Knigge 2002), or hydrodynamic numerical modeling 
(Proga, Stone & Drew 1999). A major disagreement be- 
tween the modeled line profiles and the observed ones ap- 
pears in the fit of the P-Cyg emission component. Most 
of the wind models make use the standard disk model 
(Shakura & Sunyaev 1973, Pringle 1981). Disk models 
with temperature given by the standard model arc op- 
tically thick in Balmer lines, but the outer parts could 
be optically thick in Balmer lines and optically thin in 
the continuum (Williams 1980). Following this model the 
lines must be formed in the outer parts of the disk, but the 



observations show that the disk line cmissivity arc cen- 
trally concentrated. Absorption lines arc obtained from 
optically thick disk simulations under LTE while NLTE 
effects must be considered to simulate emission lines. Ko 
et al. (1996) presented a model where the accretion disk 
vertical structure and NLTE effects are considered, where 
a chromosphere is superposed to a standard model disk. 
A selfconsistcnt model for accretion disk extended atmo- 
sphere is still an open issue. 

The main wind features are detected in UV, but wind 
effects can also appear in the optical region. In a brief sur- 
vey, Kafka & Honcycutt (2004) noticed P-Cyg profiles in 
the Hel A5876, Hel A7065 and Ha fine profiles of BZ Cam, 
Q Cyg, HR Del, DI Lac, BT Mon and AC Cnc. They pre- 
sented a few spectra of each object showing P-Cyg features 
and extended Ha line wings, that are signatures of mass 
loss by wind. Studies analyzing the behavior of these fea- 
tures with the orbital phase in the optical region were not 
presented until now. Disk spectral synthesis and eclipse 
mapping studies suggest that gas in the primary Roche 
lobe (e.g. Baptista & Catalan 2001) contribute with scat- 
tering free-free and recombination continuum photons. 

The presence of wind in CVs is a hypothesis raised to ex- 
plain the emission at low velocities seen in many Doppler 
maps and the single peaked line profiles observed even on 
systems were a double peak line profile is expected (e.g. 
high inclination systems) (Murray & Chiang 1997). 

2. Accretion Disk Simulations 

The simulation procedure is divided into two parts. 
First, a steady accretion disk emission is simulated, later, 
the flickering and/or wind is added. The simulated disk is 
geometrically thin and limited by the white dwarf at the 
inner rim and by the tidal radius at the outer rim. The 
emissivity of the disk is assumed to follow a power law 
with radial dependence 
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where Fj is each element flux, rj its radial coordinate, 9 
its angular coordinate, a the power law proportionality 
constant, b the power law index, i the orbital inclination 
and A the emission wavelength calculated from the veloc- 
ity Doppler shift calculated from equation 2. For each 
disk element the emission is calculated, the corresponding 
velocity is computed (eq. 2) and the contribution of each 
disk element is put into the spectra. 

v{r,9) = 7+ 



t^x + ^y(^-rcos(0)) +{rMd)fj (2) 
X sin(6' — 27r0)sin(i) 

In equation 2, 7 is the systemic velocity, vk the 
Keplerian disk velocity, and ax the system's stellar com- 
ponents separation. This procedure is repeated until all 
elements of the disk were computed. There is also the 
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possibility of simulate a hot spot or a boundary layer. 
These regions are geometrically restricted on the disk sur- 
face and their emissivity is given as a fraction of the total 
disk emissivity. 

Aiming to test our simulation method, a disk with an 
enhanced emission in a region similar to the hot spot at 
the outer parts of the disk and a disk with an enhanced 
emission in the inner disk parts were simulated. From the 
synthetic spectra generated by the simulations, Doppler 
maps were constructed. In the first case, the hot spot-like 
emission appears on the expected region of the synthetic 
Doppler maps, but the emission of the inner parts of the 
disk does not appear on the tomograms. From this we 
conclude that the Doppler tomography method is appro- 
priated to study emission from the outer parts of the ac- 
cretion disks but not for emission from the inner parts of 
the disk. In this later case, the information concentrated 
on a small region in position space is spread on a large 
region on the velocity space, due to the Jacobian of this 
transformation. 

2.1. Flickering Simulations 

The emission line flickering is simulated as a set of dis- 
crete flares on the accretion disk. Each flare is gener- 
ated at a random position inside a pre-defined region (the 
hot spot, for example), allowing us to simulate flickering 
from different regions of the disk or from the secondary 
star. The flickering flare has a random intensity between a 
maximum intensity (simulation parameter) and zero. This 
maximum intensity is given in disk total line flux units, 
aiming to quantify the flare energy. 

When flickering is observed in an object, the resulting 
data comprises a sum of all flickering flares events during 
each integration time. Aiming to reproduce this obser- 
vational effect on the simulations, the flickering flares are 
represented along a temporal grid with bins much smaller 
than the integration time. The sum along many bins is 
made to synthetize the effect of flickering on each spec- 
trum. As the flickering flares are not instantaneous events 
but they present a temporal evolution, the flares are sim- 
ulated as having instantaneous rise and exponential decay 

Ht) flare = I {h) flared^^ ^ (3) 

where I{t) fiare is the flare intensity on each bin, /(to) flare 
its intensity on the previous bin, r the characteristic decay 
time and t the temporal variable. The study of CVs light 
curves by Bruch (1989) indicate an assimctry on the flick- 
ering flare temporal evolution, with the flare rise faster 
than its decay. A fiducial flare description is not a critical 
issue on our simulations as the integration times are often 
longer than the individual flare timescale. 

As the flickering flares are not produced with a fixed 
periodicity, an occurrence probability is attributed to the 
flare on each time bin. This probability together with 
the integration time give us the averaged number of ex- 
pected flares on each simulated line profile. The time of 
occurrence of the flare inside a particular temporal bin is 
a random value. 

This procedure is repeated for all bins inside the inte- 



gration time. The Doppler shift of the flares are calculated 
based on the position of each flare, assuming a Keplerian 
disk, and the flare contribution is added to the line proflle. 
Aiming to simulate line profiles with noise, we introduced 
a Gaussian noise in the simulations, this noise is quantified 
by a continuum signal-to-noise ratio. The standard devi- 
ation of the noise distribution is calculated by the ratio 
between the continuum intensity and the signal-to-noise 
ratio desired for the line profile. The introduced noise is 
symmetric around zero average and it is added to each 
wavelength resolution clement of the resulting line profile. 

An averaged disk line profile and its corresponding RMS 
is shown in figure 1. Here the flickering simulations were 
held without noise. The disk was simulated on a 200x200 
bins spatial grid. The addopted simulation parameters 
were the followings: orbital period of 5 hours, primary 
mass of 1 solar mass, mass ratio of 0.7 and integrated line 
flux of 9 X 10^^^ erg s^^ cm^^. The flickering source re- 
gion is a small region restricted between 0.75 and 0.85 disk 
radii units, centered at an angle of 120? (measured from 
the line that join the stellar centres) and with an angular 
range of 1.°. This region has a location similar to the one 
of the impact of the stream on the disk considering free 
particles and a disk radius of approximately 2/3 of the 
primary Roche lobe radius. The flickering temporal evo- 
lution is calculated in a grid with 5 times more elements 
than the predicted number of flares on each integration 
time. The total decay time addopted for each flare was 1 
minute. The maximum amplitude of the flares into each 
exposition was 10% of the disk total flux. The flare oc- 
currence mean frequency is 0.025 s~^, which corresponds 
to a flare at each 40 seconds on average. The relative 
intensity between the line profile and its RMS is depen- 
dent on these parameters. The fiickering was restricted 
to a hot-spot-like region, the effect of this assumption is 
seen on figure 1, as the RMS profile lacking high velocity 
wings. The noise on the RMS profile is produced by the 
finite number of phase bins where the variance is calcu- 
lated. This behavior tends to be minimized when a large 
number of phase bins is used. 

Other simulations with flickering originating on the hot 
spot region were held. The simulation parameters were 
varied one by one to verify the effect on the synthetic 
tomograms. The addopted simulation parameters were 
the same of the last paragraph, but with an integrated 
flux of 10^^ erg s^^ cm~^, flickering maximum amplitude 
of 1% of the total line flux and including noise in the line 
profiles. 

The addopted quality criterion to quantify the detection 
of a fiickering feature was the ratio between the maximum 
intensity in the flickering occurrence region and the RMS 
of a region without flickering, affected only by the noise. 
We addopted feature detection criterion a quality 

factor [max /RMS) higher than 10. 

In the figures 2 and 3, the effect of each flickering pa- 
rameter of our simulation is presented. The error bars for 
all the simulations were calculated from the standard de- 
viation of the quality factors obtained from a series of 11 
simulations. 
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Fig. 1. Averaged simulated line profile (solid line) and its 
fliekering RMS (dotted line) line profile. The simulations pa- 
rameters are given on the text. 

From figure 2a one can see that, higher the number of 
spectra used to calculate the flickering Doppler map, bet- 
ter is the structure detection capability. The flickering 
amplitude also affects the structure detection, as features 
with higher amplitude are detected more clearly. If we 
consider that features with max/RMS> 10 are detected, 
it can be noticed that, even using a high number of spec- 
tra to calculate the maps, the flickering with lowest am- 
plitude could not be detected. On figure 3, the synthetic 
flickering maps of two particular points from figure 2 are 
shown. One can see that the feature detection is better 
for figure 36 {S/N = 20) than for figure 3a {S/N = 11). 
The flickering amplitude effect on the Doppler maps is 
shown in figure 2h. It is clear that the detection is better 
for the case of higher flickering amplitudes. One also can 
see that low amplitude flickering information is lost more 
easily than the high amplitude flickering one. The effect 
of the flare occurrence frequency is presented in flgure 2c. 
We notice a decrease on the detection quality for high fre- 
quencies, where the combination of discrete flares tends to 
form a continuum emission instead of a series of discrete 
events. It is important to observe that, in the case where 
the frequency is small than 1 flare per integration time, 
the variance calculated in each phase box may not repre- 
sent the variability properly. The integration time effect 
upon the synthetic flickering maps is shown in figure 2d. 
One can sec from this graph that the detection is better 
for small integration times (with many flares on each inte- 
gration time). For longer integration times, the flickering 
information is lost. We have not found any signiflcative 
variation of the quality factor with the decay time-scale 
of each flickering flare. 

On these simulations the flickering source can be also 
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Fig. 2. Behavior of the quality factor with (a) S/N ratio, (b) 
flickering amplitude, (c) flickering frequency and (d) integra- 
tion time. In (a) the flUed circles correspond to simulations 
with 1000 spectra and flickering amplitude 1%, the triangles 
were found with 500 spectra and amplitude of 0.5% of the to- 
tal line flux and the empty circles correspond to simulations 
with 1000 spectra and amplitude of 0.5%. 
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attributed to the reprocessing of radiation at the illumi- 
nated face of the secondary star. The application of this 
particular source model to the observations of V3885 Sgr 
was already made (Ribeiro & Diaz 2007). In this case, the 
secondary is considered spherical and the flickering source 
geometry is approximate. These assumptions are reason- 
able, due to the typical resolution of Doppler tomograms. 
The velocity field in this case is given by eq. 2 without 
the Keplerian velocity term. 

2.2. Wind Simulations 

The wind is implemented in a 3d grid, below and above 
the disk. The wind emission process is chosen to be ra- 
diative recombination, with emissivity given by 

Iwtnd -^hvaeffNj (4) 
47r 

where / is the wind emissivity, v the frequency, aeff 
the effective recombination coefficient and Ne the electron 
density. For simplicity, the wind is considered isothermal. 
As only Ha line profiles are simulated, other emission pro- 
cesses as scattering are not considered. If the temperature 
variation along the wind were considered, then a detailed 
temperature stratification on the disk must also be taken 
into account. The wind addopted geometry is biconic (fig- 
ure 4) , and its velocity law is given by the Castor & Abbott 
velocity field (Long & Knigge 2002) 

( [2GMl \ {l/Ra.r 

V,olo.dal = VO Voj (^T^Q^ (5) 

where f3 the wind terminal velocity scale factor. Mi the 
primary mass, vq the initial wind velocity, Rav the ef- 
fective acceleration scale, I is the distance of each wind 
element along the wind streamlines, a the acceleration co- 
efficient and r is the distance between a point in the wind 
and the white dwarf. The ionization was considered to be 
partial for H and He and up to the CNO elements. The 
disk angular momentum conservation through the wind is 
also considered by 

'^keplerianij'^ — '^keplerianij' d) '. 77 T (6) 

rsmya^) 

where rd is the radius of each wind element projected over 
the disk and 0^ is the elevation angular coordinate. 

As the single peaked line profiles often seen in some 
CVs could not be generated only by varying the wind pa- 
rameters, an optically thick wind was implemented. The 
line self-absorption through the wind was parameterized 
by the line optical depth t, which relates to the emitted 
radiation by 

/.(s)=/.(so)e-- (7) 

where Iv{sq) is original the wind emission and I{s) the at- 
tenuated wind emission. Scattering effects were not con- 
sidered. 

The velocity of each wind element projected on the line 
of sight is calculated. If the element is not occulted by 
the optically thick disk, its emissivity is added to the disk 




Fig. 4. Illustration of the wind biconic geometry. The white 
dwarf is the open circle at the center, the disk is the black 
thick line. The impact parameter d is the distance between 
the white dwarf and the wind focus at each side of the disk. 
The wind cone is limited by the inner and outer disk rim, and 
it is presented as the gray gradient region. The wind cone at 
the other disk side is limited by the dashed line. 

spectra. If the wind has r 7^ the absorption along the 
wind region of the disk emission is also considered. 

The result of these simulations are line profiles of 
Keplerian accretion disks with wind inside the primary's 
Roche lobe. The simulation can be repeated at many or- 
bital phases, allowing us to construct Doppler maps with 
these synthetic spectra. 

In this work the line profiles were simulated considering 
only the wind, i.e, the disk emission is not included, vary- 
ing only the optical thickness to verify the effect of this 
parameter over the line profiles (figure 5). From figure 
5 one can see that, as optical thickness is increased the 
line profile changes from double-peaked to single-peaked. 
In this simulation, we find that the line profile becomes 
single-peaked for r values higher than 10. The effect of 
optically an thick wind on the Doppler maps is presented 
on figure 6. The optically thin Doppler map (figure 6a) is 
not filled at low velocities, while the optically thick wind 
Doppler map (figure 66) presents a significant low velocity 
emission. 

The effect of the other wind simulation parameters (see 
cq. 6) on the optically thick wind line profiles are pre- 
sented on figure 7. These simulations were held with a 
fixed line optical depth t = 12.5. The impact parameter 
c?i, is what determine the collimation of the wind. Higher 
the cLtj value is, more collimated is the wind. From fig- 
ure 7a one can see that a more collimated wind produces 
a smaller flux and a narrower profile. The acceleration 
coefficient a controls the wind acceleration. In figure 76 
one can see that a more accelerated wind generates most 
intense line profiles. One interesting point about this pa- 
rameter is that, for low values of a, a deficit in the emission 
at the red side of the line profile is perceived. This can be 
understood as the wind region opposite to the observer 
(redshifted mainly by the poloidal velocity component) 
having its radiation more attenuated by the absorption 
than the side in the direction of the observer (blueshifted) . 
If only the Keplerian wind velocity is considered, the at- 
tenuation would be equal in both red and blue sides. The 
effect of terminal velocity scale factor /3 on the line profiles 
(figure 7c) is that a higher value of j3 yields line profiles 
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Fig. 3. Two synthetic flickering Dopplor maps for 1000 spectra, flickering maximum amplitude of 0,5% the total line flux and (a) 
S/N = 11 and (6) S/N = 20. Notice that the feature detection is much better on the second case. The crosses are, from top to 
bottom, the secondary's center of mass, the LI Lagrange point, the system's center of mass and the primary's center of mass. 
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Fig. 5. Effect of optical depth on the wind line profiles. For 
optical depths higher than approximately 10 the lines are sin- 
gle peaked. 

with higher flux. A higher line flux is also the result of 
an increase of the effective acceleration scale Rav (figure 
7d) and a decrease on the initial wind velocity vq (figure 
7e). In these last three cases the line profile shape does 
not present a significant change with the variation of each 
parameter, only the line flux does. 

Another physical scenario that produces single-peaked 
line profiles and consequently low velocity emission on 
Doppler maps is the emission by stationary material in- 
side the Roche lobe. In this case, the only velocity com- 
ponent is orbital and the density of the material inside the 



Roche lobe is assumed to be constant. The line profiles 
generated by emission of stationary material inside the 
primary's Roche lobe are similar to the ones produced by 
high optical depth winds. In the present simulations wind 
and/or emission by stationary material can be considered. 

A test application was performed to verify the plau- 
sibility of the wind hypothesis to produce low velocity 
emission on observed Doppler maps. Wind Doppler maps 
were simulated with some different values of line optical 
depth. From each map a velocity profile was obtained by 
calculating the mode along rings centered on the primary. 
These velocity profiles were compared to the V3885 Sgr 
Doppler map velocity profile (Ribeiro & Diaz 2007) (fig- 
ure 8). From figure 8 one can see that the V3885 Sgr Ha 
velocity profile has a shape similar to the simulated wind 
map velocity profile for high line optical depth The veloc- 
ity of the peak does not coincide because it depends on 
the dynamical parameters of the wind. 

We also analyzed the effect of the line optical depth 
on the radial emissivity profiles. We generated a series of 
simulated wind line profiles calculated with some fixed val- 
ues for the line optical depth. Using those scries Doppler 
maps were computed for each case. The disk contribution 
to the total emissivity was not included. The Doppler 
maps were transformed to the position space assuming a 
Keplerian velocity regime to obtain the radial emissivity 
profiles. The modal emissivity value was then calculated 
along several rings centered on the primary. Finally, a 
power law was fitted to these points (figure 9). The in- 
dexes of these radial power laws as a function of the line 
optical depth are given on table 1. 

From table 1 we can conclude that the radial emissivity 
gets less centrally concentrated as wc have higher line op- 
tical depths in the wind. From figure 9 one can see that, 
as expected, the flux is smaller for higher optical depths, 
but it is rather low for the r/i„e = 15 case. This is expected 
as the emission measure gets too small in the high optical 




Table 1. Power law indexes of the radial emissivity 
wind profiles as a function of the line optical depth 
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Fig. 7. Effect of the wind parameters on tiie wind line pro- 
files, for a line optical depth r = 12.5. The parameters are 
(a) impact parameter, (6) acceleration coefficient, (c) termi- 
nal velocity scale factor, (d) effective acceleration scale and 
(e) initial wind velocity. The high frequency features that 
appear mainly in the line wings arc effect of the finite grid 
addoptcd in the simulations. 
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Fig. 8. Velocity profiles from (a) Hq {solid curve) and Hel 
A6678 (dashed curve) V3885 Sgr Doppler maps and from (6) 
simulated wind Doppler maps. On (b) the profiles are pre- 
sented for some line optical depth values: r = {solid black 
curve), T = 10 (dotted black curve), t = 15 (dashed black 
curve), T = 20 (solid gray curve) and t = 25 (dashed gray 
curve). The V3885 Sgr Ha line velocity profile has a shape 
similar to those simulated with 15 ^ r "J 20. 
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Fig. 9. Radial emissivity profiles from simulated Ha wind 
Doppler maps. The data are displayed with alternated colors 
for each t value for clarity 
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depth regime. The flux is also dependent on the other pa- 
rameters considered on the simulations, as the electronie 
density or the dynamical parameters of the wind. The 
next step of our work is analyze the behavior of the sim- 
ulated velocity and radial emissivity profiles with these 
parameters and compare the real data profiles with the 
simulated ones to constrain these dynamical parameters 
using a multiparametric optimization. 

3. Discussion 

In all flickering tomography studies the integration time 
must be small enough to preserve the flickering informa- 
tion. On the other hand, it must be long enough to ensure 
a good S/N ratio. One limitation of the Doppler tomogra- 
phy technique is that the flickering from the inner parts of 
the accretion disk could not be mapped with this method, 
as the information from a proportionally small accretion 
disk region is spread over a large region on the Doppler 
map. In addition, most of the emission from the inner 
parts of the accretion disk falls on the line wings, the pro- 
file region that is most affected by noise and continuum 
subtraction errors. 

One possible extension of the flickering simulation study 
is to compare the observed flickering maps and simulated 
ones aiming to constrain the flare energy. The difficult of 
this study is that the flickering simulation has many pa- 
rameters and that the noise is superposed to the intrinsic 
stochastic variation on the observed data. The flickering 
in the present work is simulated considering only one oc- 
currence frequency, more frequencies can be included in 
future studies. 

When the flickering from the secondary star is simu- 
lated, it is supposed to be due to the flickering from the 
disk being reprocessed on the illuminated face of the com- 
panion. To produce a more fiducial simulation, the effect 
of the illumination itself must be considered on the code. 
The effects of the illumination over the disk line profile is 
not included in the present version of the code. 

On both the wind and flickering simulations, the disk 
emission is given by a power law emissivity radial profile. 
The present simulation code can be combined with an disk 
model code to improve the description of the disk emission 
properties. One obstacle to do this is that longer compu- 
tational times will be necessary to simulate the spectra 
used on the construction of the synthetic Doppler maps. 
On the wind simulations, the wind was considered isother- 
mal. In the case of a disk model combined with the code, 
the temperature stratification of the wind can be other- 
wise included. 

The SW Sex stars are cataclysmic variables with pe- 
culiar characteristics. Among others, they present single 
peaked line profiles where a double peaked line profile is 
expected and their Doppler maps are dominated by emis- 
sion at low velocity. This is occasionally observed even on 
systems which are not classified as SW Sex CVs. Many 
scenarios have been proposed to explain this behavior: ac- 
cretion disk winds (Horne 1986), magnetic accretion and 
mass loss via a propeller mechanism (Horne 1999), a gas 
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stream that does not impact the disk completely at the hot 
spot region and flyby over the disk (Hellier & Robinson 
1994), accretion on the magnetic poles of a synchronous 
primary (Casares et al. 1996) and shock along the stream 
trajectory over the disk (Groot, Rutten & van Paradijs 
2000) . The present simulations and their comparison with 
observed Doppler maps may help to identify the origin of 
SW Sex phenomenon. In particular, the reconstructions 
of a disk with an optically thick wind suggest that disk 
mass loss should be at least part of the solution to the SW 
Sex problem. 

4. Conclusions 

A method for simulating accretion disk Doppler maps 
including wind and flickering is presented. The detec- 
tion quality of the flickering raising from a small region of 
the disk was analyzed. High S/N and high time resolu- 
tion spectra are needed to obtain information from flicker- 
ing tomograms. The information from low amplitude and 
high frequency flickering is lost first. The flickering infor- 
mation is also lost if long integration times are used. The 
flickering from the inner part of the accretion disk could 
not be mapped with the Doppler tomography technique. 

Wind simulations were also held, and the effect of each 
simulation parameter over the wind line profiles is pre- 
sented. It was necessary to consider an optically thick 
wind, with optical depth higher than 10 to generate single 
peaked Ha line profiles. Diffuse emission inside the pri- 
mary Roche lobe is other mechanism that generate these 
single peaked profiles and centrally concentrated Doppler 
maps. 
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